Measuring Quark-Gluon-Plasma Thermalization Time with Dileptons 
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We calculate the medium dilepton yield from a quark-gluon plasma which has a time-dependent 
momentum-space anisotropy. A phenomenological model for the hard momentum scale, Phard(T), 
and plasma anisotropy parameter, £(t), is constructed which interpolates between longitudinal free 
streaming at early times (r <C n sa ) and ideal hydrodynamic at late times (r 2> n so ). We show 
that high-energy dilepton production is sensitive to the plasma isotropization time, Ti so , and can 
therefore be used to experimentally determine the time of onset for hydrodynamic expansion of a 
quark-gluon plasma and the magnitude of expected early-time momentum-space anisotropics. 
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An important question facing experimentalists and 
theorists working on heavy-ion experiments ongoing at 
the Relativistic Heavy Ion Collider (RHIC) and planned 
at the Large Hadron Collider (LHC) is at what time is 
it justified to assume that the matter created can be de- 
scribed using simple hydrodynamics. At RHIC energies 
it has been found that for pt < 2 GeV the elliptic flow of 
the matter created is described well by models which as- 
sume ideal hydrodynamic behavior starting at very early 
times t < 1 fm/c [1-4]. Since then further refinements 
including viscous corrections have become available and 
indications are that the viscosity is constrained to be 
small [5-7] and results remain consistent with early ther- 
malization of the quark-gluon plasma (QGP). However, 
given the complexity of solving three-dimensional viscous 
hydrodynamic equations coupled to a late-time hadronic 
cascade [8, 9] it would be nice to have an independent 
way to determine the time at which a QGP begins to 
undergo hydrodynamic expansion. In this letter we pro- 
pose to use high-energy dilepton yields as a function of 
both pair invariant mass and transverse momentum to 
experimentally determine QGP isotropization time. 

The use of ideal hydrodynamics to describe matter re- 
quires at a minimum that the matter be isotropic in mo- 
mentum space [10]. In practical applications it is also 
necessary to impose an additional constraint (equation- 
of-state, conformal invariance, etc.) in order to close the 
resulting system of hydrodynamic equations. When an 
equation-of-state is applied it is implicitly assumed that 
the system is isotropic and thermal so that t; so = T t herm- 
For simplicity we will also identify these two time scales. 
Estimates from pcrturbative QCD for the thermalization 
time of a QGP range from 2 — 3 fm/c [11-13]. Recently it 
has been shown that plasma isotropization is accelerated 
by unstable plasma modes induced by the rapid longitu- 
dinal expansion of the QGP fireball [13, 14]; however, it 
is still not known by precisely how much. 



In most phenomenological treatments of the QGP it 
is assumed that the plasma thermalizes rapidly with 
T iso = Ttherm on the order of the parton formation time. 
However, given the rapid longitudinal expansion of the 
matter this seems like a rather strong assumption and 
one would like to know the effect of momentum-space 
anisotropics on experimental observables. 

Absent a precise dynamical picture of the first few fm/c 
of the QGP's lifetime we propose a simple phenomeno- 
logical model for the time-dependence of the plasma 
momentum-space anisotropy, £ = \{Pt)I{Pl) — 1> an d 
hard momentum scale, phard- We then use this model to 
explore the effect of early-time plasma momentum-space 
anisotropics on high-energy dilepton production. To ac- 
complish this we introduce two parameters: (1) t iso which 
is the proper time at which the system begins behav- 
ing hydrodynamically and (2) 7 which sets the sharpness 
of the transition to hydrodynamic behavior. For times 
greater than the parton formation time, To, but short 
compared to T{ so we will assume that the system is longi- 
tudinally free streaming and for times long compared to 
Tiso that it is expanding hydrodynamically. To estimate 
the parton formation time we use the nuclear saturation 
scale, r ~ Qj 1 [15]. For RHIC energies Q s ~ 1.5 GeV 
and for LHC energies Q s ~ 2 GeV implying that To ~ 0.2 
fm/c and To ~ 0.1 fm/c, respectively. 

For our final results we present the invariant mass 
and transverse momentum dependence of medium dilep- 
ton production. We compare our results with estimates 
of other relevant sources and background processes and 
demonstrate that increasing the QGP isotropization time 
from T iso = t to Tj so — 2 fm/c significantly increases 
medium high-energy dilepton production. We present 
detailed calculations of this enhancement and show that 
at LHC energy it leads to an experimentally measurable 
effect on dilepton production between 3 GeV < M, pr < 
8 GeV. 
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DILEPTON RATE 

The leading order medium dilepton production rate 
comes from the annihilation process qq — ► . Medium 
dilepton production is computed by folding this rate 
together with the expected local quark and anti-quark 
phase space distribution functions [16, 17]. In this work 
we allow these phase space distribution functions to 
be anisotropic in momentum-space and then model the 
time-dependence of the anisotropy. Assuming azimuthal 
symmetry about the beam (z) axis the time-dependence 
can be parameterized as 

/{<?,<?} ( k ' T ) = f{q,q}-}so (\A 2 +£(t) k 2 z , PhardO)) , (1) 

where phard(T) is a time-dependent hard momentum scale 
and £(t) is a time-dependent parameter reflecting the 
strength of the local momentum-space anisotropy [18]. 
We will further assume that /i SO is a Fermi-Dirac distri- 
bution and that f q = fq. 

Details of the analytic and numerical calculation of 
the ^-dependent rate will be presented elsewhere. Note 
that for isotropic systems there are appreciable next-to- 
leading (NLO) order corrections for E/T < 1 [19-22]. 
When we are sensitive to areas of phase space where there 
are large NLO corrections we will apply A-factors to our 
estimates as indicated. 

SPACE-TIME MODEL 

We now construct a model which interpolates between 
early-time longitudinal free streaming and late-time ideal 
hydrodynamic expansion. Given (1) the medium parton 
energy density can be factorized as 

£(Phard,f) = J 7-^3 Pf(P,0 = £o(Phard) K(Q , 

(2) 

where K{£) = [l/(£ + 1) + arctanVf/V?] /2 and 
£o(Phard) is the energy density resulting from integra- 



tion of the isotropic quark and anti-quark distribution 
functions appearing in Eq. (1). 

Longitudinal Free Streaming Limit: For a longitudi- 
nally free streaming plasma £fs( t ) — ( T / r o) 2 — 1 and 
Phard is constant and equal to the initial average hard 
momentum scale in the plasma; therefore, £fs{ t ) — 
£(Phard> £fs(t))- Assuming an isotropic plasma at r = To 
this results mphard(r) = T and lim T>To £ F s -> £o (to/ t ) 
where To is the initial plasma "temperature" . Note that 
assuming a formation time of To = 0.1 fm/c at r = 1 
fm/c we have £fs — 100. 

Ideal Hydrodynamic Expansion: For a plasma which is 
undergoing ideal longitudinal hydrodynamic expansion 
we have £(r) =0 by assumption. Additionally, since the 
system is thermal we can identify the hard momentum 
scale with the plasma temperature so that Phard(T") = 

1/3 

T(t) = To (to/t) . Correspondingly, we have £hydro = 

£ (Vt) 4 / 3 . 

Interpolating Model: In order to construct a model 
which interpolates between longitudinal free streaming 
and hydrodynamic expansion we introduce a smeared 
step function A(r) = (tanh [j(t — T iso )/T ] + 1) /2. This 
allows us to model the time-dependence of £ and phard as 

£{t) = £ FS (r) [U(r)/U(T )} i/3 , 
Phard(r) = T [U(t)/U(t )} 1/3 , 

H(t) = a 2 ^)) _ i , ( 3 ) 

where U{r) = [^(a? so - l)] 3A(r)/4 (a iso /a) A(T) , a = t/t 
and ajso = Ti so /To. The power of TZ in hi keeps the energy 
density continuous at r = t- iso for all 7. 

When t <C t- 1so we have A — > and the system is 
longitudinally free streaming. When r ^> r iBO then A — > 1 
and the system is expanding hydrodynamically. In the 
limit 7—5-00, A — > Q(t — Ti so ). In Fig. 1 we plot the 
time-dependence of £, phard, and £ assuming 7 = 2 for 
four different plasma isotropization times corresponding 
to a iso e {1,4,8,16}. 
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RESULTS 

To obtain the final expected dilepton yields we inte- 
grate the annihilation rate over t G { T o,i~f} and rj £ 
{—2.5,2.5} with parameters specified by Eq. (3) and r/ 
set by Phaid^Tf) — T c . In this letter we will present ex- 
pected e + e~ yields resulting from a Pb-Pb collision at 
LHC full beam energy, yfs = 5.5 TeV. At RHIC energies 
sensitivity to r iso is reduced due to the poor signal-to- 
background ratio for medium dileptons. Predictions for 
Au-Au at RHIC energy will be presented elsewhere. 

In order to facilitate comparison with previous works 
we take r = 0.088 fm/c, T = 845 MeV, T c = 160 
MeV, and Rt — 7.1 fm [22]. Here we assume that when 
the system reaches T c all medium emission stops. The 
addition of mixed and hadronic phase emission is not 
included in the present work since the kinematic range 
we consider is not sensitive to the late-time evolution of 
the system (see Fig. 5). In Figs. 2 and 3 we show our final 
predicted e + e~ yields as a function of invariant mass and 
transverse momentum along with predicted yields from 
other sources. Predictions for Drell-Yan, heavy quark, 
and jet conversion yields were supplied by the authors of 
Ref. [22]. 

As can be seen from Fig. 2 there is a significant vari- 
ation of the medium dilepton yield when varying the as- 
sumed plasma isotropization time from 0.088 fm/c to 
2 fm/c. When an isotropization time of 2 fm/c is as- 
sumed we see that medium dileptons become as impor- 
tant as Drell-Yan and jet conversion. The reason for the 
enhanced production is that longitudinal free streaming 
preserves more transverse momentum than an hydrody- 
namically expanding plasma. However, all three contri- 
butions are down by an order of magnitude from the 
expected background coming from semileptonic heavy 
quark decay. In practice this would mean that exper- 
imentalists would have to subtract this background to 
10%. We note that an indefinitely longitudinally free 
streaming plasma (r; so — > oo) produces less low energy 
(M,pt < 2 GeV) dileptons due to the rapidly dropping 
parton densities as £ — > oo. 

As we show in Fig. 3 as a function of pt the medium 
contribution dominates the expected Drell-Yan and jet 
conversion sources for all pr < 6 GeV. If an isotropization 
time of 2 fm/c is assumed then the medium dileptons 
dominate out to px ~ 9 GeV. This dominance means 
that it should be possible to use dilepton production to 
determine much-needed information about quark-gluon 
plasma initial conditions at LHC. As can be seen from 
Fig. 3 at pr — 5 GeV the expected medium dilepton 
yield varies by nearly an order of magnitude depending 
on the assumed plasma isotropization time. This level of 
variation will hopefully be measurable at LHC. 

To illustrate the dependence on the model parameter 7 
in Fig. 4 we have plotted medium dilepton yields obtained 
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FIG. 2: Dilepton yields as a function of invariant mass with 
a cut pt > 8 GeV. For medium dileptons we use 7 = 2 and 
ri so is taken to be either 0.088 fm/c or 2 fm/c. A iVfactor of 
1.5 was applied to account for NLO corrections. 



i-9 
o 



> 
O 



Q. 

T3 

■o 



— H~ 


Drell Yan 




Jet Conversion 




Medium (r is0 = 0.088 fm/c) 




Medium (r ]s0 = 2 fm/c) 




6 

Pt [GeV] 



10 



FIG. 3: Dilepton yields as a function of transverse momentum 
with a cut 0.5 < M < 1 GeV. For medium dileptons we use 
7 = 2 and Ti so is taken to be either 0.088 fm/c or 2 fm/c. A 
A"-factor of 6 was applied to account for NLO corrections. 



assuming t- iso = 0.5 fm/c and Ti so = 2 fm/c. The central 
values obtained are with 7 = 2 and the error bars come 
from variation of 7 in the range 0.05 < 7 < 10. As can 
be seen from this Figure between 3 and 8 GeV there is 
little sensitivity to the parameter 7. 

In Fig. 5 we show snapshots of the fraction of medium 
dileptons produced as a function of transverse momen- 
tum. What this figure shows is that by 4 fm/c yields 
in this kinematic regime are saturated. At 1 fm/c ap- 
proximately 94% of all pt — 4 GeV dileptons have al- 
ready been produced as well as 98% of the pt = 5 
GeV pairs. This highlights the sensitivity of this observ- 
able to early-times after a heavy-ion collision and jus- 
tifies neglecting the effect of transverse expansion and 
mixed/hadronic phases when considering dileptons in 
this kinematic regime. 
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FIG. 4: Dilepton yields as a function of transverse mo- 
mentum. Shown are yields obtained assuming Ti so = 0.5 
fm/c and 2 fm/c with error bars indicating model variation, 
0.05 < 7 < 10. Cuts and Ti'-factor are the same as in Fig. 3. 
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FIG. 5: Fraction of dileptons produced at r £ {1, 2, 4, 16} 
fm/c assuming Ti so = 0.5 fm/c. Cuts are the same as in Fig. 3. 



CONCLUSIONS AND DISCUSSION 

Based on Figs. 3 and 4 it should be possible to measure 
Tj SO at LHC energies using dilepton production in the 
kinematic range 3 < pr < 8 GeV. Additionally via our 
model for £(t) given in Eq. 3 determining Ti so provides 
an estimate of the maximum amount of momentum-space 
anisotropy achieved during the lifetime of the QGP. 

The effect of varying r iso is also large in the dilepton 
spectra vs invariant mass as shown in Fig. 2, however, 
Drell-Yan and jet conversion production can be up to 10 
times larger than medium production making it difficult 
to measure a clean medium dilepton signal. 

Our chief uncertainty is the NLO order corrections to 
dilepton production incorporating anisotropics. These 
corrections are particularly important for low- mass dilep- 
ton production. At LHC energy it is possible to re- 
duce sensitivity to these NLO corrections by placing cuts 
M, pt > 2 GeV. Another uncertainty comes from our as- 
sumption of chemical equilibrium. Naively finite chemi- 



cal potentials should affect isotropic and anisotropic plas- 
mas equally so one expects that although the total yields 
could change one would still see a sensitivity to the as- 
sumed isotropization/thermalization time. At leading or- 
der in the quark fugacity, X q , the ratio of the isotropic to 
anisotropic rates should be independent of X q [23] . 

Future work will study the effect of finite quark chemi- 
cal potentials, collisional broadening of the parton dis- 
tributions, and the possibility of late-time persistent 
anisotropies (finite viscosity). In addition, models such 
as (3) can be used to assess the impact of momentum- 
space anisotropies on other observables. 
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